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1. Introduction 

The lack of a systematic calculational scheme applicable to the infrared sector of 
Quantum Chromodynamics (QCD) has motivated the advent of models aspiring to 
provide phenomenologically viable links between asymptotic freedom and confinement, 
by enriching perturbation theory with judiciously selected nonperturbative information. 
An important source of such information is provided by the dispersion relations. The 
latter, being based on the "first principles" of the theory, furnish the definite analytic 
properties with respect to a given kinematic variable of a physical quantity under 
consideration. The basic idea behind the so-called "analytic approach" to Quantum 
Field Theory (QFT) is to supplement perturbation theory, and in particular the 
renormalization group (RG) formalism, with the nonperturbative information encoded 
in the relevant dispersion relations. Specifically, the perturbative solutions of the RG 
equations possess unphysical singularities in the infrared domain, a fact that contradicts 
the general principles of local QFT. The analytization procedure amounts to the 
restoration of the correct analytic properties for a physical quantity at hand, by forcing it 
to satisfy the Kallen-Lehmann spectral representation (see equation (jl))). This method 
was first proposed in the framework of Quantum Electrodynamics (QED) and applied 
to the study of the invariant charge of the theory pQ. Later on, these considerations 
were generalized to the case of the non-Abelian theories, and the analytic approach 
to QCD [2] emerged. This approach has been successfully applied to the study of the 
strong running coupling [210], perturbative series for the QCD observables [I], and some 
intrinsically nonperturbative aspects of the hadron dynamics [HI |3] . Some of the main 
advantages of the analytic approach to QCD are the absence of unphysical singularities 
(by construction), and a fairly good higher-loop and scheme stability of the results 
obtained. 

A central quantity within the aforementioned approach is a(k 2 ), the running (or 
"effective") coupling (charge) of QCD. Clearly, this quantity can be reliably computed 
only in the ultraviolet region, where perturbation theory is applicable, and must be 
modelled in the infrared domain, where perturbative methods break down, and one 
eventually encounters the unphysical singularities, such as the Landau pole. The 
analytic approach to QCD eliminates such artefacts and provides a concrete analytic 
expression for the running coupling in the infrared, maintaining at the same time 
the standard asymptotic behaviour in the ultraviolet. However, since the analyticity 
requirement can be incorporated into the RG formalism in various ways, two main 
pictures have emerged. Specifically, if the analyticity condition is imposed directly on 
the perturbative running coupling a s {k 2 ) one arrives at the model of [2] with finite 
universal infrared limiting value given by ct(0) = Ait / (3q. If, instead, the analyticity 
requirement is imposed on the corresponding (3 function, one obtains the running 
coupling of 3J which is singular ("enhanced") at k 2 = 0. The latter invariant charge 
has proved to be able to describe a number of the strong interaction processes both of 
perturbative and of nonperturbative nature [3J 13] . Undoubtedly, it would be interesting 
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to further study the physical implications and phenomenological possibilities offered by 
the infrared enhanced analytic invariant charge. The primary objective of this paper is 
to examine the compatibility of the infrared enhanced analytic running coupling with 
chiral symmetry breaking (CSB), and its impact on the solutions of the Schwinger- 
Dyson (SD) equation controlling the dynamical generation of quark masses. 

The way the running coupling enters into the standard SD equation for the quark 
propagator ( "gap equation" ) is rather well-known. Since QCD is not a fixed point theory, 
the QED-inspired gap equation must be modified, in order to incorporate the asymptotic 
freedom. The usual way of accomplishing this eventually boils down to the replacement 
1/k 2 — > a(k 2 )/k 2 in the corresponding kernel of the gap equation. The inclusion of a(k 2 ) 
is essential for arriving at an integral equation for the quark propagator S(p) which is 
well-behaved in the ultraviolet. Indeed, the additional logarithm in the denominator 
of the kernel due to the strong running coupling a(k 2 ) improves the convergence of 
the integral. However, since the perturbative form of a(k 2 ) diverges as l/ln(/c 2 /A 2 ) 
when k 2 — > A 2 , a physically motivated model for the infrared behaviour of a(k 2 ) is 
needed. The infrared enhanced analytic charge represents a good candidate for such a 
purpose, since, as has been explained in detail in [SI E] , it combines asymptotic freedom 
and confinement behaviour in a single expression. This is to be contrasted with the 
standard method of introducing asymptotic freedom and confinement (see, e.g., 6j), 
whereby the running coupling is composed by two separate, and theoretically rather 
uncorrelated, contributions (see further discussion in Section HJ). 

The paper is organized as follows. In Section El we briefly review the most salient 
features of the analytic approach, with particular emphasis on the running coupling 
displaying infrared enhancement. In Section El we go over the usual assumptions and 
approximations entering into the derivation of the standard gap equation. In Section 0] 
we study the asymptotic behaviour of the dynamical quark mass function, by converting 
the integral equation into a differential form. In Sectional the integral equation is solved 
numerically, and the results are presented. It turns out that, with four active quarks, 
the experimental value of the pion decay constant f n may be obtained if the QCD mass- 
scale A, which is the only free parameter available, acquires the rather elevated value of 
about 880 Me V. Finally, in Section |H1 we discuss the results and present our conclusions. 

2. Analytic invariant charge in QCD 

In this Section we give a brief summary of one of the models for the QCD analytic 
invariant charge [3]. This model shares all the advantages of the analytic approach 
and it was successful in the description of a wide range of QCD phenomena both of 
perturbative and of intrinsically nonperturbative nature 0. Furthermore, it is of a 
particular interest to note that this model has recently been re-derived [7j, proceeding 
from completely different motivations. 

As has already been commented in the Introduction, the basic idea behind the 
analytic approach to QFT [H Ej is to supplement perturbation theory with the 
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nonperturbative information encoded in the relevant dispersion relations. The latter 
are based on "first principles" of the theory, and provide one with the definite analytic 
properties of a given physical quantity with respect to a proper kinematic variable. 
In practice the "analytization procedure" [2] amounts to the restoration of the correct 
analytic properties for a given quantity F(k 2 ) by imposing the Kallen-Lehmann integral 
representation 1 

{F( e)} =f4U. (i) 

L Jan Jq o -f k 2 

Here the spectral function g(a) is determined by the initial (perturbative) expression 
q(o) = lim \F(-a - is) - F(-a + is)] , (2) 

with a > 0. 

The model for the analytic invariant charge 0E1 that we will study in this paper is 
obtained by imposing the requirement of analyticity (PJ) on the perturbative expansion 
of the RG (3 function 

j'+i' 



d lnajg(/i 2 
d In jj? 



47T 



(3) 



In this equation a^(/i 2 ) is the £-loop analytic invariant charge, a^(/x 2 ) denotes the 
£-loop perturbative running coupling, (3j stands for the (3 function expansion coefficient 
(Po = 11 — 2n { /3, /3i = 102 — 38n f /3, ... ), and ra f is the number of active quarks. It 
is worth noting that prescription (J3J) differs from that of the original Shirkov-Solovtsov 
model [2] , where the requirement was imposed directly on the perturbative running 
coupling a s {k 2 ) (discussion of this issue can also be found in [HI El E])- 
At the one-loop level the RG equation @ can be solved explicitly [Hj: 

The solution to equation (jHJ) can also be represented in the form of the Kallen-Lehmann 
integral 

mm-*?. 



Po Jo cr + z 
where the one-loop spectral density is 

P (1) (-) = (l + ^)pr^, (6) 

V cry In a + n z 

and the explicit expression for the £-loop (cr) can be found in (HI |H] . 

The invariant charge flU) possesses a number of appealing features. First of all, 
it has the correct analytic properties in the k 2 variable, demanded in equation ((TJ, 
namely, it has the only cut k 2 < along the negative semiaxis of real k 2 . Then, it 

1 A metric with signature (—1, 1, 1, 1) is used, so that positive k 2 corresponds to a spacelike momentum 
transfer. 
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Figure 1. The analytic invariant charge a4n(fc 2 ) (see equation Q) and the 
perturbative running coupling a^ik 2 ) at the one-loop level (solid and dashed curves, 
respectively). Here rif = 4 active quarks are assumed. 



contains no adjustable parameters 2 . Thus, similarly to the perturbative approach, the 
QCD scale parameter A remains the basic characterizing quantity of the theory. In 
addition, the invariant charge (@J) incorporates the ultraviolet asymptotic freedom with 
the infrared enhancement in a single expression (see Figure [Q), that plays an essential 
role in applications of the model in hand to the description of the quenched lattice 
simulation data (see 0E1 for the details). It is worth mentioning here that the singular 
behaviour of the strong running coupling a(k 2 ) when k 2 — > is also supported by a 
number of studies of the SD equations (see, e.g., papers [111111111121113] and references 
therein). Moreover, the invariant charge (jSJ) displays a good higher loop and scheme 
stability, and it has proved to describe a number of strong interaction processes (e.g., 
confining static quark-antiquark potential, inclusive r lepton decay) in a self-consistent 
way. The detailed analysis of the properties of the analytic running coupling © and its 
applications can be found in [31|5J|B]. 

Given the characteristic features of the analytic invariant charge mentioned above, 
it would be interesting to study its influence on phenomena particularly sensitive to 
the low-energy dynamics. To that end, in the following three sections we study how 
the analytic invariant charge (0) affects the the mechanism of CSB and dynamical mass 
generation for the quarks, through the study of the SD equations governing the dynamics 
of the quark propagator [T^l HS1 HZl EHl HH1 1201 [2H 1221 - 



2 It should be noted that the Shirkov-Solovtsov running coupling [2] has no adjustable parameters, 
either. So, both these models are the "minimal" ones in this sense. 
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3. The gap equation 

Throughout this Section, we will work exclusively in Euclidean space. According to the 
usual conventions, the starting point is to express the fully dressed quark propagator in 
the following general form [T8] : 

S-\p) =tf + m + E(p) = i$A(p 2 ) + B(p 2 ), (7) 

where m is the bare quark mass and E(p) is the quark self-energy. 

Before proceeding, it is worthwhile to make some additional comments on the above 
functions. Since we will be concentrated in the case where one does not have explicit 
CSB, i.e., the bare mass mo = 0, the quark mass is generated only through dynamical 
effects. In this case, CSB takes place when the self-energy S(p) develops a nonzero 
value. Alternatively, one may define the quark mass function M(p 2 ) in terms of the 
functions A(p 2 ) and B(p 2 ), as M(p 2 ) = B(p 2 )/A(p 2 ); then CSB is considered to occur 
when B(p 2 ) ^ 0. 




E S V 



Figure 2. The SD equation JHJ for the quark self-energy. The black blobs represent the 
fully dressed quark and gluon propagators and the white one is the proper quark-gluon 
vertex. 

The quark SD equation is represented schematically in Figure El and can be written 

as 

E(p) = ~/j ^J»S(q)T„(q,p)D^(k), (8) 

where we have used that J2 a A a A a = 4/3, A a being the Gell-Mann matrices, and k = p—q. 
According to this equation, the self-energy is dynamically determined in terms of 
itself, the full gluon propagator, denoted by D ab u {k) = 5 ab D^(k), and the full quark- 
gluon vertex T u (q, p). Of course, both D°%,(k) and T u (q,p) obey their own complicated 
SD equations, a fact which eventually makes the use of simplifications unavoidable and 
further modelling of the unknown functions involved. 

A common approximation employed in the literature (see, e.g., CHI E3J is to 
neglect the ghost contributions in the quark SD equation, whose effects are supposed to 
be partially accounted for by the fully dressed gluon propagator and the full quark- 
gluon vertex. This assumption leads to a theory with Abelian-like characteristics, 
where the usual non-Abelian Slavnov-Taylor identities are replaced by QED-like Ward 
identities [23] • I n particular, for the quark-gluon vertex we have 

tk»T,(p,q) = S- 1 (p)-S- 1 (q). (9) 



7 



The imposition of this so-called "Abelian approximation" gives rise to further 
simplifications. Specifically, due to the (assumed) validity of equation (jUJ), the usual 
QED identity Z\ = Z2, where Z\ and Z2 are the renormalization constants for the 
fermion-boson vertex and fermion wave function, respectively, is restored. This, in turn, 
allows the definition of a RG invariant quantity, to be denoted by a(k 2 ), which is the 
exact analogue of the QED effective charge, namely 



g 2 D, u {k) 



k^ky 



A-Ka(k 2 



(10) 



The exact transversality of the right hand-side of equation (|TU|) is the result of working 
in the Landau gauge. A different choice of gauge would have resulted in an additional 
longitudinal term of the form ^k^k^/k 4 , where £ is the usual gauge-fixing parameter 
(£ = corresponds to the Landau gauge, £ = 1 to the Feynman gauge). Note that 
the nature of this extra term is purely tree- level, i.e., there is no higher-order dressing 
involved. 

Furthermore, the Ward-identity © motivates the use of the time-honored "Gauge 
Technique" [2*lj . Specifically, a nonperturbative Ansatz for the vertex r^p, q) in terms 
of S(p) is postulated, based on the requirement that it should satisfy, by construction, 
equation (JHJ). Evidently, such a construction leaves the transverse part of T^p, q) 
unspecified; the usual argument around this ambiguity is that, in a theory with a mass- 
gap, the transverse parts of the vertex are sub-leading in the infrared, and have little 
or no consequence on CSB (see, e.g., [21] )• In the rest of our analysis we will use for 
^ flip j l) the simple Ansatz proposed in jTUl EH] 

iT M (p, q) = iA(p 2 )^ + ^{1 [A(p 2 ) - A(q 2 

Choosing the Landau gauge leads to the further simplification 

g 2 D^(k)T u (q,p) = g 2 D^(k)A(q 2 ) lu , 

since, in that case, the gluon propagator is completely transverse. 

Substituting equations (J7J) and (fT2*|l into quark gap equation (jSJ) 
commonly used coupled system for the quark self-energy [TOJ HH1 I2H] 



6 + [B(p 2 )-B(q 2 )}}. (11) 



(12) 



one arrives at the 



A(p 2 

and 
B(p 2 ) 



V 



d 4 q A7ia{k 2 



(2 



7T 



k 2 



p-q + 2 



p-k q-k s 



A\q 2 ) 



k 2 I q 2 A 2 (q 2 ) + B 2 (q 2 



d 4 q A7ia(k 2 ) A{q 2 )B(q 2 



(13) 



(14) 



(2tt) 4 k 2 q 2 A 2 (q 2 ) + B 2 (q 2 ) 
Note that the angular integration can be easily evaluated in the last two equations, by 
resorting to the usual angle approximation for the running coupling [HH 12*3*1 127)] 



a 



((P - <?) 2 ) « 0(p 2 - q 2 )a(p 2 ) 



+ 



p 2 )a{q 2 



(15) 



In particular, with such an approximation, the angular integral for the Dirac-vector 
component of the quark self-energy vanishes, leading automatically to A{p 2 ) = 1. 
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Therefore, we can straightforwardly relate B(p 2 ) to the dynamical mass M(p 2 ). Then 
the coupled system formed by equations (fT5|l and (JHJ), reduces to one single equation 
namely [2H] 

~a(x) r x yM(y) f°° a(y)M(y) 



M{x) = - 

7T 



Jo v + M 2 (v) V Jx v + M 2 (v) V 



(16) 



x Jo y + M 2 (y) J x y + M 2 (y) 
where M{x) = M(xA 2 )/A, x = p 2 /A 2 , and y = q 2 /A 2 . 

4. Asymptotical behaviour of the mass function 

The main effect of implementing the substitution described by equation (jlOjl at the level 
of the gap equation is to transfer our ignorance regarding the behaviour of the gluon 
propagator into the nonperturbative structure of the invariant charge a(k 2 ). This latter 
quantity can be modelled more directly, essentially because it enters more naturally 
than the gluon propagator in the parametrization of the low-energy QCD data. In the 
following sections, we will study in detail the solutions obtained from the quark gap 
equation (fTT)|) after using as a(k 2 ) the one-loop analytic running coupling a^(fc 2 ) (jlj). 

The inclusion or not of appropriately modelled confinement effects has an important 
impact on the type of solutions that one obtains from the gap equation. In fact, it has 
often be claimed in the literature that, if such effects are not included, one may not 
encounter non-trivial solutions to the gap equation at all, i.e., CSB does not occur 
(see, e.g., [27]). The usual way of accounting for confinement effects at the level of 
the gap equation is to insert a gluon propagator of the form A 2 /k 4 , the (spacial) 
Fourier transform of linearly rising quark-antiquark potential [2H]- Evidently, this 
expression fails to capture asymptotic freedom in the ultraviolet domain, whose effects 
are separately supplemented through the inclusion of the corresponding perturbative 
contributions. Thus, the "standard" way of describing both effects is through a running 
coupling of the form 

/.2\ C Air 1 k 2 

a — (fc) = 7 + ^R7T7)' Z = W (17) 

where C and r are dimensionless constants (see, e.g., |TU] ) . Usually, C is treated 
as an adjustable parameter, to be determined in such a way as to reproduce the 
correct phenomenological values for the quark masses, pion decay constant, and chiral 
condensates. On the other hand, r plays the role of an infrared regulator; it is often 
treated as an arbitrary parameter, but in a more complete, physically motivated picture 
of QCD, it has been identified as a dynamically generated gluon "mass" |2Z1 EH EHj- 
Note that, if C = 0, the logarithmic term on the right hand-side of (|17|) must overcome 
comfortably an infrared critical coupling of about 1.2, in order to obtain from the gap 
equations phenomenologically interesting solutions. This may or may not be possible, 
depending on whether r is considered as a free parameter, or if some physical arguments 
constrain its possible range of values, as is the case with the gluon "mass" . 

A definite advantage of the analytic invariant charge (k 2 ) (Jlj), compared to (JT7J), 
is the simultaneous incorporation of asymptotic freedom and infrared enhancement into 
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a single expression, of concrete theoretical origin, namely the analyticity properties of the 
theory. In particular, note that, unlike a coupling of the type given in ([17)1. the running 
coupling (@J) contains no adjustable parameters. This theoretically appealing feature, 
is, of course, much more restrictive when one handles phenomenologically relevant 
quantities. It is worth noting also that the enhancement displayed by the analytic 
coupling a^(k 2 ) has been shown to correspond to the confining static quark-antiquark 
potential with a quasilinear rising behaviour at large distances. Namely, V(r) ~ r/lnr, 
when r — > oo, with r being the dimensionless distance between quark and antiquark, 
see E] for the details. 

Before proceeding with the numerical solution for the SD equation (fTKj) for all range 
of momenta, it would be interesting to gain some explicit insight into its behaviour in 
the deep infrared region. On general grounds, given that the running coupling (jH) 
diverges at the origin, one might also expect a similar behaviour from the solutions 
of equation (fTTIj) . In fact, this is what happens when (fTtjj) is solved using the running 
coupling (fTTj) . which is also singular in the infrared limit [TUj . 

In order to obtain the low-energy asymptotic behaviour of the dynamical mass 
function A4(x), it is convenient to cast the integral equation (jlfij) into a differential 
form, by differentiating both sides twice with respect to x (see also [TU| 123]) 



d 

dx 



d_ 

dx 



xM(x) 
a(x) 



d_ 

dx 



1 a(x)M(x) 

7T X + M 2 (x) 



a(x) 

By making use of the explicit form of the one-loop analytic invariant charge (jlj) one can 
reduce this equation in the limit x — > to 

M(v A^ + 2 i) M(v) = h v ^-°°' (19) 

where v = Inx. In deriving (j!9|) we have neglected the factor z in the numerator 
of (@J), since in the limit considered it is subleading next to 1. The above equation 
can be solved analytically through successive iterations, giving rise, as expected, to 
divergent solutions, which are formally expressed in terms of powers of v and lnt>. The 
first iteration, corresponding to the leading infrared behaviour of the dynamical mass 
function A4(x), is obtained by omitting the second-order derivative in equation (|T9*|) . 
Then, the solution is 



M{x) ~ J— In | lnx |, x^O. (20) 
V Po 

The divergence of the dynamical quark mass function Ai(x) in the low-energy domain 
x — > can be interpreted as a hint for confinement, see also discussion in [TTH 131 j . 

It is worthwhile to emphasize again that the obtained behaviour ({20]) is restricted 
to the deep infrared domain (i.e., for values p 2 A 2 ). Subleading corrections to this 
solution may be systematically obtained from equation however, they are of the 
same order as the terms neglected when arriving from (|18|) to (|19|). and are therefore 
of little usefulness. It is also interesting to compare equation (|2T)|) to the corresponding 
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solution obtained by making use of the running coupling (|17|) . The leading low-energy 
behaviour in that case is instead Ai KTW (x) ~ y2C\lax\/7T, when x — > 0. Thus, we 
infer that for the case of the analytic invariant charge the infrared singularity of 
the dynamical mass function M.(x) is much milder than in the case of the running 
coupling (fTTj) . 

For the sake of completeness, we finish this Section by reporting the ultraviolet 
asymptotic expression for the dynamical quark mass Ai(x). Since we are considering 
the case of exact chiral symmetry (i.e., no bare mass, m = 0), the conservation of the 
axial-vector current eventually leads, for sufficient large momenta, to 

Mix) ~ — (lnx) A -\ x^oo, (21) 
x 

where A = 12/(33 — 2ni) is the anomalous dimension of the mass, and D is a 
constant independent of the renormalization point and directly related to the quark 
condensate pffi 



5. Numerical solution for the mass function 



The numerical solution for the dynamical quark mass function Ai(x), obtained solving 
directly the integral equation (jTBj) , is presented in Figure El Indeed, one can see 
numerically a soft increase of the mass function Ai(x) when x — > 0, as suggested by the 
leading behaviour (|2Uj) extracted from the differential equation analysis. Regarding the 
ultraviolet region, we note that choosing sufficiently large values for the ultraviolet cutoff, 
the results obtained are independent of the latter, since the integrals in equation (fTfij) 
are ultraviolet convergent. Typically, our solutions are evaluated within a momenta 
window of sixteen orders of magnitude (10~ 8 < x < 10 8 ), which is sufficient to ensure 
their stability. 

Now, with the solution for the dynamical mass at hand, we can relate the value of 
the QCD scale parameter A to the pion decay constant f n , defined as the axial- vector 
transition amplitude for an on-shell pion. This can be accomplished by making use of the 
method developed by Pagels, Stokar j^l] , and Cornwall jM] , which is a generalization of 
the Goldberger-Treiman relation when the momentum carried by the pion is different 
from zero: 



J 7T 



3 A 2 f°° yM(y) 



4vr 2 Jo [y + M 2 (y)Y 



M(y)- 



ydM{y) 



dy. (22) 



2 dy 

Thus, one is able to fix the value for A, the only adjustable parameter in this 
analysis, by requiring that the pion decay constant should assume its measured value of 
f n = 93 MeV [33] . For the more favorable case of n? = 4 active quarks this procedure 
results in the estimate A = 880 MeV. 

The inclusion of the higher loop corrections to the analytic running coupling a an (k 2 ) 
(see equation (jHJ) is not expected to alter the qualitative picture obtained above. This 
is so because the most intrinsic feature of this charge, namely the infrared enhancement, 
persists after the incorporation of the loop corrections; however, the type of singularity 
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Figure 3. The dimcnsionless quark dynamical mass M{x) — M(xA 2 )/ A, obtained 
from equation (|l(jfl using the one-loop analytic invariant charge Q for Uf = 4 quark 
flavours. 



displayed at the origin becomes slightly milder than in the one-loop case. Therefore, 
we anticipate that the confining behaviour of M.(x) will also persist, but with a weaker 
infrared singularity then that of equation fl20(l . On the other hand, in this case one 
would expect a higher estimate for A, given that the loop corrections lower the value 
of a an (k 2 ) in the entire range of momenta [Sj. This fact, in turn, will lead to smaller 
values for Ai(x), making the saturation of the right hand-side of ()22j) more difficult, 
and forcing A to assume even higher values. 

We finish this Section by commenting on the veracity of the angular approximation 
given in equation (fT5|). and the dependence of the obtained results on it. It would 
be certainly of interest to establish whether the divergent nature of the solutions 
persists, or is an artefact of the aforementioned approximation. Indeed, one could 
envisage the possibility that the simultaneous solution of the system (fl"3|) and (JHJ) 
might actually lead to a finite expression for M(x), as x — > 0, despite the fact that the 
kernel is divergent at the origin, due to the enhanced form of the running charge (J3J) 
employed. In order to address this point in some detail, we have not resorted to the 
approximation of equation (fl3|) . but have instead performed the angular integration 
numerically, and subsequently attempted to solve the resulting coupled system. The 
numerical integration of the final (momentum) integrals requires the introduction of an 
infrared regulator [36J; we have regulated the kernels by carrying out the replacement 
a(k 2 )/k 2 —>■ a(k 2 )/(k 2 + fi 2 ). Evidently, the solutions for the dynamical mass function 
depend now explicitly on the regulator /i, and one should study the behaviour of 
A4(x, fi 2 ) in the limit fi 2 — > 0. Our numerical analysis revealed that there is a 
critical value of /i 2 , of about jj 2 ~ 1CT 7 A 2 , below which no convergent solutions to the 
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system (|13|) and (|14j) may be found. Although no inescapable conclusion may be drawn 
from this fact, we interpret this breakdown as a strong indication that the resulting 
solutions diverge as fi 2 — > 0, in qualitative agreement with what was found when the 
angular approximation of equation (fTB]) was used. 

6. Conclusions 

In this paper we have studied the compatibility between the infrared enhanced analytic 
charge of QCD and chiral symmetry breaking, through the detailed analysis of a 
standard form of the gap equation for the quark propagator, where this former coupling 
was incorporated. It turned out that, due to the infrared enhancement of the running 
coupling employed, the solutions found for the dynamical quark mass function M(p 2 ) 
were infrared divergent. Following standard arguments we have interpreted this 
divergent behaviour as an indication of confinement. The final upshot of this analysis 
was that the aforementioned analytic charge is able to break the chiral symmetry, 
furnishing a reasonable ratio between the QCD scale parameter A and the pion decay 
constant 

To be sure, the value of 880 MeV obtained for A (with four active quark flavours) 
appears elevated when compared to the "standard" value of A of about 350 MeV quoted 
in the literature (see [US] and references therein), but also when compared to the values 
for A obtained within the analytic approach itself by resorting to different methods [Sj, 
ranging between 500-600 MeV. 3 We emphasize, however, that the main purpose of the 
analysis presented is not so much to extract an accurate value for A, but rather check to 
what extend two a priori different methods, the analytic approach and the SD equations, 
may coexist in a complementary and qualitatively consistent picture. In that sense 
we consider the outcome of the present work encouraging, especially when taking into 
account the theoretical uncertainties intrinsic to both methods, and the fact that, unlike 
the majority of existing models, in our case A is the only adjustable parameter available. 

It would certainly be worthwhile attempting to improve the above picture by 
incorporating into the spectral density defining the analytic charge (see equation (jH))) 
contributions from nonperturbative effects (e.g., the operator product expansion (see 
also [HB]) an d the nonlocal chiral quark model [US])- We hope to be able to report 
progress in this direction in the near future. 
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